PROCEEDINGS 

OF SCIENCE 



in 
o 
o 

(N 

Or 
<D 

m 

(N 



On 
O 

in 
o 



X 



UV-filtered overlap fermions 



Stephan Durr 

Bern University 

E-mail: |durr (AT) itp . unibe . ch 

Christian Hoelbling 

Wuppertal University 
E-mail: |hoelbling (AT) physik . uni-wuppertal . de 



Urs Wenger 

NIC/DESY Zeuthen 



E-mail: urs . wenger (AT) desy . de 



We discuss the kernel spectrum, locality properties and the axial-vector renormalization constant 
of UV-filtered overlap fermions. We find that UV-filtered overlap fermions have a better condi- 
tioned kernel, better locality and an axial-vector renormalization constant closer to 1 than their 
unfiltered counterparts, even if the shift parameter p is simply set to 1 . 
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Figure 1: Degree of the Chebychev polynomial required to achieve a minimax accuracy of 10~ 8 on a 16 4 
lattice at j3 = 6.0. With p = 1 a single HYP step brings for the massless overlap operator a speedup factor 
~ 4, but a fair comparison is to the situation with p = 1 .4 without filtering, and then the factor is ~ 2. 



1. Introduction 

Overlap fermions ^ have many theoretically desirable features, but they are computationally 
rather demanding. There are suggestions in the literature on how to reduce the computational cost 
by using a more elaborate kernel together with fattened gauge links [§, ||, Q]. More recently, it 
has been observed that UV-filtering alone has a similar effect [§, ^, |7[ |9[ and maintains the 
0{a 2 ) Symanzik class. In this note we investigate the effectiveness of UV-filtering in the quenched 
approximation. We use the Wilson gauge action and employ APE [11] or HYP [ |l2| ] smearing. For 
details of the implementation and simulation parameters we refer the reader to []9| |lOj ]. 



2. Kernel spectrum 

The overlap operator is constructed from the hermitean Wilson operator //w,m = 75Dw.n1 via 



D ov = p (1 + sign(# w _ p )) 



(2.1) 



where p is the shift parameter. In a practical implementation one has to find a method to approx- 
imate the sign function in ( |2.1[ ) over the entire eigenvalue spectrum of Hy/-p- Therefore, it is 
essential for an efficient implementation that the kernel operator //w,-p is well conditioned. 
In the free theory the condition number C(p) of |fiw,-p I i s given by 



C(P) 



±f- for < p < 1 
§Ejf for 1 < p < 2 



(2.2) 



which has a minimum value 7 at p = 1. In the interacting theory the upper edge of the spectrum at 
8 — p is only slightly affected for typical p\ On the other hand the absence of a lower bound on the 
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Figure 2: Condition number 1/e = C(p) of \Hy/.- p \ versus p for j3 = 5.84 (left) and j3 = 6.0 (right). 



spectrum can result in a dramatic increase of the condition number C(p). Since there are standard 
techniques to treat the lowest few of these eigenmodes exactly, the relevant issue is the onset of the 
bulk of these modes, and it turns out that the latter gets significantly lifted by UV-filtering 

In order to translate the condition number into computational cost, we choose to present the 



degree of the Chebychev polynomial [13] necessary to reach a given precision in the minimax norm 
after treating the 14 lowest eigenmodes exactly. As one can see from Fig. [T], the number of forward 
applications of Hyj needed to obtain a given precision is significantly reduced by UV-filtering. 

More details are shown in Fig.^. A rather conservative filtering recipe like applying 1 step of 
APE link fattening seems to render the kernel far more benevolent; the condition number remains 
comfortably small for a large range of shift parameters, including the canonical choice p = 1 . 



3. Locality 

The "spread" inherent in the gauge link fattening raises the question whether a UV-filtered 
overlap operator is less local than the original thin link version. It has already been observed in [^] 
that the UV-filtering actually proves beneficial to the locality properties of the overlap operator. 

Fig. H shows the localization of the overlap operator at /3 = 6.0. Following the convention of 



[ 14 ] we plot the maximum over the 2-norm of D ov rj at x with tj a normalized 8 -peak source vector 



at the point in the lattice versus the "taxi driver" distance d\ = \ \x\\\ to the location of the 5-peak 



/(^)=SUp(j|(A>vT?)(*)|b 



d x . (3.1) 



Looking first at the unfiltered operators (black/dark diamonds and crosses) one finds the well known 
result that (at this j8) adjusting p to a value around 1.4 lets f{d\) fall off steeper than keeping 
the untuned value p = 1. The interesting observation is that a single HYP step together with 
p = 1 (red/light squares) results in an even steeper descent than the unfiltered version with p = 1.4 
(which was chosen to nearly optimize the locality of the unfiltered operator). The last curve shown 
(red/light pluses) indicates that one should not attempt to combine the filtering with a p value that 
would be optimal for the unfiltered operator. 
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Figure 3: Localization of the overlap operator at j3 = 6.0 without filtering and after 1 HYP step, for p = 1 
and p = 1 .4. A single HYP step proves more efficient than optimizing p . Filtering and p > 1 should not be 
combined; for the filtered operator the untuned choice p = 1 is reasonable (but still not optimal). 



4. Renormalization 

One interesting aspect of UV-filtered fermions is their improved renormalization behavior. It 



has been realized some time ago [ ]15| , [16J that in general one-loop corrections to renormalization 
constants are substantially smaller in the case of UV-filtered fermions. 

To investigate this behavior in a nonperturbative context we choose to determine the axial- 
vector renormalization constant Za from the AWI, following the method of 0, [l8|]. Our results 
are collected in Tab. |l| [gj. It is interesting to note that for the untuned p = 1 thin link operator Za is 
large even at /3 = 6. With p = 1 .4 and the same coupling Za is much closer to 1, indicating that p > 
1 is not only essential for the locality of the unfiltered D ov but also beneficial to its renormalization 
constants. Note, however, that even the 1-fold APE smeared operator with the untuned p = 1 has 
an axial- vector renormalization constant which is closer to 1. When going to the coarser j8 = 5.84 
lattice Za of the thin link operator with p = 1 increases drastically. At this point it becomes essential 





5.66 5.76 5.84 6.00 


Z^ one ,p = 1.4 [19] 

z none p = L g [2Q] 

Z] APE ,p = l 
Z] APE ,p = l 

z a hyp ,p = i 

Z] HYP ,p = l 


ill-def. ill-def. 7.06(73) 3.145(94) 
— — ~ 1.71 1.553(02) 
_ _ ~1.44 — 
2.57(7) 1.90(2) 1.66(02) 1.452(04) 
1.55(3) 1.33(1) 1.23(01) 1.160(06) 
1.44(2) 1.28(1) 1.22(01) 1.153(03) 
1.21(1) 1.13(1) 1.10(01) 1.072(02) 



Table 1: Za with several filterings. Note that [|19| g(J] give Za values at slightly different /3 only. The precise 



numbers are Z™ ne Q3 =5. 8458, p = 1.4) = 1.710(5) [jig and Z™ ne (j3 = 5.85, p = 1.6) = 1.443(5) [|2 
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Figure 4: Axial-vector renormalization constant Za versus j3 . The thin link data at p = 1 .4 are from [11 



filtering level 


none, p = 1.4 [[TJ] 1 APE 3 APE 1 HYP 3 HYP 


C\ +c 3 

C3 


1.66(16) 0.36(32) 0.03(20) 0.15(15) 0.03(13) 
5.54(4) 5.53(2) 5.51(4) 5.49(4) 5.49(5) 



Table 2: The 1-loop coefficient c\ + cj, and pole location C3 of a fit (4. 1) for different levels of UV-filtering 
Note, that the analytic value for c\ + C3 for thin link overlap at p = 1 .4 is c\ + C3 = 0.722408 [EU E2n. 



for the thin link operator to properly tune p . Going to even coarser lattices we were unable to obtain 
a signal for Z A with the unfiltered operator and p = 1 (for details see [^]). 

In Fig. |] the Za with different filterings are plotted vs. j8 . The lines represent fits of the form 



Z A (x) 



1 +C1X + C2X 
1 — C3X 



(4.1) 



with x = 1//3 and all fits have a reasonable % 2 . To compare to the thin link p = 1.4 case we use 
data from Q19[]. In principle, these curves contain two pieces of information. The asymptotic slope 
c\ +C3 for x —s- predicts the perturbative 1-loop coefficient for Za and the pole in ( fOl ) [i.e. the 



value of C3] predicts the coupling where the ansatz ( |4.1| ) breaks down. With the current quality of 
our data all asymptotic slopes - except for the thin link operator - are compatible with zero and we 
are unable to make a quantitative statement about the perturbative 1-loop coefficients (see Tab. |2[). 
Still, on a qualitative level the filtered one-loop coefficients are smaller than both the analytic result 
and the fit to nonperturbative data for thin link overlap with p = 1.4. This suggest a much better 
behaved perturbative series for the filtered overlap operator [ph. For the location of the pole the 
results in Tab. ^ indicate that this quantity is hardly changed by the filtering. This suggests that the 
(practical) range of validity of the perturbative expansion is barely affected by the filtering 
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5. Conclusions 

We have presented some evidence that UV-filtered overlap fermions might bring substantial 
technical advantages over their thin link counterparts. In particular, a forward application of the 
massless overlap operator requires fewer applications of the hermitean Wilson operator, the locality 
of the overlap operator is generally improved with UV-filtering and the axial- vector renormalization 
constant is much closer to 1. In addition, there is no need to tune the kernel shift parameter p. The 
canonical choice p = 1 is satisfactory even on rather coarse lattices. 
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